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Abstract: It is common practice in passive microwave remote sensing (microwave 

radiometry) to express observables as temperatures and in the case of polarimetric 
radiometry to use what are called “Modified Stokes Parameters in Brightness 
Temperature” to describe the scene. However, definitions with slightly different 
normalization (with and without division by bandwidth) have appeared in the literature. 
The purpose of this manuscript is to present an analysis to clarify the meaning of terms in 
the definition and resolve the question of the proper normalization. 

I. Introduction 

It is common practice in passive microwave remote sensing (microwave radiometry) to 
express observables as temperatures, a practice apparently inherited from radio 
astronomy [ la, 2a]. In the case of polarimetric radiometry it is common to use the 
notation of Stokes parameters but with the correlations normalized to units of temperature 
(e.g. [3]). The terminology has recently been clarified in the “Standard Terminology for 
Microwave Radiometry”, a monograph of definitions prepared by the radiometry 
community, where the normalized correlations are called “Modified Stokes Parameters in 
Brightness Temperature” [4, 5], The definition in Chapter III of this monograph is: 


T v 

[X 2 /qkBw] < | E v 1 2 > 

(la) 

T h 

[k 2 /qkB w | < | E h | 2 > 

(lb) 

t 3 

[T 2 /r|kB\vl 2 Re <E V E h *> 

(lc) 

t 4 

[7 2 /qkB w ] 2 Im <E V E h *> 

(Id) 


where X is the wavelength, k is the Boltzmann’s constant, B w is bandwidth and q is the 
wave impedance of the medium. As will be pointed out below, there is ambiguity with 



the meaning of expressions such as < | E v 1 2 >. There also appears to be disagreement 
regarding the form of Equations 1 and versions have appeared in the literature with and 
without the presence of bandwidth, , B w , on the right hand side [6, 7, 8], The purpose of 
this manuscript is to present an analysis, in the context of radiometry, to clarify the 
meaning of the quantities on the right-hand side of Equations 1, Our goal is to resolve 
the question of the proper normalization in these expressions and to add insight into the 
definitions presented in the “Standard Terminology for Microwave Radiometry” [4]. 

One approach for obtaining Equations 1 is to model the radiometer response along the 
lines of classical radio astronomy (i.e. in terms of the incident brightness) and then use 
arguments from electromagnetic theory to relate the brightness to the electromagnetic 
fields. This two-step approach is outlined in Appendices A and B, respectively. A 
slightly different approach is adopted here in an effort to relate the results more directly 
to polarimetric measurements as used in contemporary microwave remote sensing 
radiometry. The analysis begins in Section II with the output of an idealized correlation 
radiometer. Section III rewrites the output as equivalent temperature and Section IV puts 
the results in the form of Equations l. An interpretation in terms of familiar results is 
given in Section V. We conclude that when < | E p 1 2 >/q is equivalent to the source 
brightness at polarization “p”, then the bandwidth, Bw , should not appear in Equations 1. 

II. Radiometer Response 

Consider the idealized “correlation” radiometer shown in Figure 1. It consists of two 
paths whose signals are multiplied together and averaged. The multiplication is 
“coherent” (i.e. phase and amplitude are measured). For the purpose of this investigation, 
the details of the implementation with real hardware (e.g. [9]) is omitted and all the 
operations are assumed ideal: The antenna is lossless and its output voltage, e p (t), is the 
input to an ideal receiver (no noise, no cross-coupling, ideal operations). In this idealized 
case, h p (t) denotes the effective impulse response of the arm with polarization, p = (p,q), 
and the last step after the multiplication (operationally a low pass fdter) is assumed to be 
sufficiently narrowband to be represented by an integral. An idealized receiver is 
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appropriate here because the goal is to determine the form which the input takes after 
processing rather than the limits that real hardware place on estimating this quantity. 

A correlation receiver has been taken as the starting point to be general and also to 
explicitly include the processing involved in polarimetric radiometry. A total power 
radiometer is a special case when the two paths are identical. In the more general case, 
the system represents a single antenna with two different polarization ports (e.g. p = h 
and q = v) as is appropriate for polarimetric radiometry. In principle, the schematic could 
also represent measurements with two different antennas as would be the case, for 
example, in an interferometer or synthetic aperture radiometer [10]. The analysis in this 
case follows as outlined below, but one must explicitly include the relative distance 
(spacing) between the antennas [11]. Including this term and realistic constraints [11,12] 
complicates the analysis and since it is not pertinent to the issue at hand, this case is not 
included here 






Fig. 1: Receiver architecture. The output signal, e p (t), from each antenna 
is the input to an ideal correlating radiometer receiver. The antenna is 
described by its complex voltage patterns and the receiver path is described 
by its impulse response, h p (t). 


The output of the receiver in Figure 1 can be written: 

T 

V P q(t) = (1/2'T) J [e P (t-t’)*h p (l-t’)] [e q (t-t')*h q (t-t’)] dt’ (2) 

- T 
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Equation 2 is the result of three operations: a) The response of each arm (p,q) of the 
receiver which is given by the convolution e(t)*h(t); b) the product of the signal out of 
each arm which is the product of the two convolution pairs (i.e. the terms in square 
brackets); and c) the time average of the product. The last step is represented by the 
integral over the parameter t’. The parameter, T, represents the equivalent time constant 
of the last stage of the receiver in Figure 1 and the integration is over a window of width 
2T centered at t. In Equation 2, the asterisk indicates a convolution but in all the 
work to follow it will indicate the complex conjugate. Replacing the terms in the 
integrand of Equation 2 with their Fourier transform and using the fact that they are real 
quantities, one obtains: 

T 

V P q(t) = (l/2-2)J JdvJdv’ E p (i>)E q *(u’) Hp(\))Hq*(i)’) exp[-j27t(u-u’)(t-t’)] dt’ (3) 

- T 

Ep(u) and H p (n) denote the Fourier transform of the time domain functions, e p (t) and h p (t), 
respectively. (See Appendix B for a definition of the Fourier transform pair.) Following 
convention, E p (o), can be written in terms of the incident electric fields, 2p(Q, v). and the 
complex “voltage” patterns, J%p(Q, u), of the antennas [13]: 

Ep(n) = VA P | £p(Q, v)j? p (fl, n) dO (4) 

4k 

where A p is the effective area of the antenna [lb, 2b] and Q is solid angle: D = [0,(p]and 
dfi = sin(0)d0d(p. Substituting Equation 4 into Equation 3, one obtains: 

T 

V P q(t) = (1/22) J { jj V(ApA q ) ! J -£p(Q, v) v') JfyQ, «) ... 

- T 

... H p (n)H q *(n , )exp[-j27r(n-n’)(t-t , )l d£2 d£2’ do dn’ }dt’ (5) 
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Since the application is to remote sensing of natural sources, the incident field will be 
modeled as a random process. Consequently, the radiometer output will also be a random 
process. The goal is to compute the expected value of the radiometer output, < V pq (t)>: 

T 

<V pq (t)> = (1/27) J{JJ V(ApA q ) j J ... 

-T 

... H p (n)H q :!: (n’)exp[(-j23i(i)-u , )(t-t , )]dQ dQ’ do du’ } dt’ (6) 

To simplify the analysis, it is assumed that the process is stationary and that the fields, 
'£p(Q, o), arriving at the antennas from different directions are uncorrelated, common 
assumptions in remote sensing of “natural” scenes [14]. These assumptions can be 
represented formally by writing: 

< £p(Q,u) V(Q’^ , )> = S pq (n,u) 8(v-v’)5(Q -Q ’) (7) 

where S pq (n,n) is the cross-spectral density associated with the incident electromagnetic 
field. (See Appendix C for a discussion of the meaning of S pq (f2,o) and the use of Fourier 
transforms in this context.) Using Equation 7, one obtains: 

<V pq > = V(ApAq) J J Spq(Q,u) ;3p(Q, t>) j3q*(Q,u) H p (o)H q *(u) di2 dv (8) 

Because of the assumed stationarity, the expected value of the radiometer output, 

<V pq (t)>, does not depend on time, and the notation has been simplified to reflect this. 

When the passband of the receiver, H Piq (o), is narrow compared to the scale of change 
with frequency of the spectral density and the antenna, one has approximately: 

(9a) 

(9b) 


V(A p A q ) J S pq (an 0 )^p(n,no)^*(auo)dQ I H p (u)H q *(u) dv 
B w H° pq ]V(A p A q ) j S pq (n,i)o) ^p(Q, u 0 ) u 0 ) d^ } 
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where t»o is the center frequency of the radiometer receiver and 


B w = { J H p (d) H q *(u) dn } / H p (u 0 ) H q *(u 0 ) (10) 

H° pq = Hp(i)o) H q *(n 0 ) (11) 

The bandwidth, Bw , defined in Equation 10 is the common definition of “noise- 
equivalent bandwidth” [2c, 4]. It is implicit in the definitions above that the gain of the 
system is a maximum at the center frequency, i>o. 

III. Antenna Temperature 

The convention used in radiometry is to express the output as a temperature, T, by 
equating the output power to the noise power available from a resistor at this temperature. 
At microwave frequencies (i.e. in the Rayleigh- Jeans approximation) one has: [15,16] : 

W = kTBw- In order to obtain the radiometer output in units of power, notice that 
<V pq >/q has units of power. The choice of the impedance, q = Ve/p, for the 
normalization is justified by the application to the observation of electromagnetic fields 
and Poynting’s theorem (e.g. see Appendix B). Dividing both sides of Equation 9 by q 
one has: 

<V pq >/q = kTBw (12) 

= V(A p A q )Bw / [S pq (H)/qJ ^,(n)^*(fl)dn (13) 

For convenience, in Equation 13 the system gain has been set to unity (H° pq = 1) and the 
explicit dependence of S pq (Q) and A p , q (Q) on frequency Do has been omitted. Equation 
13 is the sum of the incident power over all directions weighted by the effect of the 
antenna. From Equation 13, it can be seen that the dimensions of S pq (£2)/q are 
Watts/(meter 2 -steradian-Hz). These are the dimensions of “brightness” [la]. (In optics 
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this quantity is called “spectral radiance” and in microwave radiometry the name 
“spectral brightness” has also been used [2d] 

To obtain results in conventional form, the effective area in Equation 13 is replaced by 
the antenna beam “solid angle” [lb]: 

A p = A, 2 /Q B p (14) 

where X is evaluated at the center frequency no and Qb p is the integral over the 
normalized antenna pattern. Equating the right-hand sides of Equations 12 and 13 and 
using Equation 14, one obtains: 

T A pq = 1 1/ V( £2b p O eq 3 I 1 (A 2 /pk) S pq (Q) jfyQ) dfi (15) 

The term, T A pq , in Equation 15 is a measure of the output power of the radiometer and is 
called the “antenna temperature”. 

IV. Modified Stokes Parameters 

In the special case p = q. Equation 15 becomes: 

T a pp = (l/n B p) J (X 2 /pk) Spp(I2) Pnptfl) da (16) 

where P np = J4p(Q, u) 34p*(f2,u) is the normalized antenna power pattern [2b]. Comparing 
Equation 16 with Equation A6 (Appendix A) one sees that the term (AT/qk) S PP (Q) in 
Equation 16 plays the role of the brightness temperature of the source: 

Tpp(O) = brightness temperature of the source at polarization, p 

= (A 2 /r|k) Spp(H) (17) 
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Alternatively, one could use the discussion below Equation 13 or in Appendix B 
identifying S PP (Q) )/r| with brightness. Then, assuming an unpolarized source and using 
the Rayleigh- Jeans approximation for applications to observations in the microwave 
spectrum, one has: 

S pp (f2) /r| = Brightness at polarization, p 

= kT pp (Q)/X 2 (18) 

which is the same as Equation 17. The factor of 2 in the Rayleigh- Jeans approximation is 
not present above because this is the power at a single polarization. 

The definition of brightness temperature can be generalized by adapting the relationship 
in Equation 17 to the general case of mixed polarization: 

Tpq(Q) = (l 2 /X]k) Spq(Q) (19) 

The temperature on the left hand side of Equation 19 is a brightness temperature (a 
characteristic of the source) in contrast to Equations 15 or 16 where the temperature on 
the left hand side represents the output of the measurement system (antenna temperature). 

Finally, adopting the notation of modified Stokes parameters for brightness temperature 
as used in polarimetric radiometry [3], one can write Equations 19 in the form of a four- 
element vector: 


T pp 


T v 


S pp (£2) 

T qq 

= 

T h 

= (A. 2 /rjk) 

Sqq(£2) 

2Re T p q 


t 3 


2 Re S P q(Q) 

21m T pq 


3 


„ 2 1m S pq (Qj_ 


( 20 ) 
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The notation [T v , Th T3 T4] in the middle set of brackets in Equation 20 is what is called 
in radiometry, “Modified Stokes Parameters in Brightness Temperature” and assumes 
that v = p and q = h. 

V. Generalized Brightness 


It is also possible to use the relationships above (e.g. Equations 18 and 19) as the basis 
for defining a generalized brightness temperature for sources with mixed polarization: 
B pq = Sp q (n) /q. Following the procedure above, one can define a 4-vector brightness: 



Bp? 


“ s pp (n, d) “ 


Bqq 


S qq (Q, 0 ) 

= 


= 1 / T] 



2 Re B pq 


2 Re S pq (fi, u) 


2 Im B pq 


2 Im S pq (G, t>) 


*» MM 


M MM 


( 21 ) 


The results can also be written in matrix form: 


B pp 

CD 

X! 

1 


S pp (H, D) / T| 

S p q(Q, 0) / q 

_ Bpq* 

Bqq 


Spq(f 2 ,u)* / q 

S qq (Q, U) / q 


(22) 


The matrix defined in Equation 22 has the form of the coherency matrix defined in optics 
[17a]. The terms in Equation 22 have dimensions of brightness. However, normalizing 
by the factor B 0 = ( Vs pp S qq )/ q yields a matrix of dimensionless elements which is 
identical to the “degree of coherence” defined in optics [17a]: 

[4 pq = [B pq J/B 0 (23) 
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Following the optics definition, the special case of unpolarized radiation (i.e. natural, 
thermal radiation), is defined to be the case when |J pq = 8 pq (i.e. zero when p 7 s q and 1 
when p = q). From thermodynamic arguments in this case (Planck law and the Rayleigh- 
Jeans approximation) one has that: B 0 = kT/). 2 . Assuming unpolarized radiation, and 
using Equation 23 and Equation 22, one has: 

Bpp + Bqq = [ S PP (Q, u) + Sqq(Q, u) ] / q (24a) 

= 2 k T /X 2 (24b) 

The term on the right in Equation 24b is just the Rayleigh-Jeans approximation for the 
brightness of an unpolarized source. From the discussion above, it is clear that it applies 
to unpolarized, “natural” radiation and represents the total power, the sum of both 
degrees of freedom (both polarizations). 

The same information is contained in the modified Stokes parameters. Using Equations 
20 and 22, it is also possible to express the brightness in terms of the modified Stokes 
parameters: 


IBvh] 


(k / X 2 ) 


T v 

(T 3 - i T 4 )/2 


(T 3 + i T 4 )/2 
T h 


(25) 


Comparing the expression above with case of plane waves [17a] suggests some important 
special cases: 

1. Natural radiation: T 3 = T 4 = 0 and T v = Th; 

2. Linearly polarized radiation: T 4 = 0; 

3. Reflected natural radiation: T 3 = T 4 = 0andT v # Th. 

Item 3 above is the circumstance encountered when microwave radiation from a thermal 
source such as the Sun or Moon is reflected from a surface such as the Earth to the 
observer. This is a common occurrence in microwave radiometry and points out the need 
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for a more generalized form for scene brightness than the classical Rayleigh-Jeans 
approximation: B = 2kT/I . It is possible in this case that T 3 ^ 0 and/or that T 4 f 0. 

The former can occur due to polarization rotation such as occurs in the ionosphere (e.g. 
Faraday rotation) or due to simple misalignment of the antenna polarization vectors. The 
later (T4 f 0) is possible with scattering from rough surfaces when polarization mixing 
occurs [18].” 

VI. Discussion 

Equations 20 are the “Modified Stokes Parameters in Brightness Temperature”. Writing 
them in the format of Equations 1 , one obtains: 


T v 

[X 2 /tfk] 

<S VV (£2,u)> 

(26a) 

Th = 

[X 2 /r]k] 

<Shh (£2,u)> 

(26b) 

t 3 

[X 2 /t]k] 

2 Re <S V h(£2,u)> 

(26c) 

t 4 

[X 2 /rjk] 

2 Im <S V h(£2,t>)> 

(26d) 


The dependence on £1 and u on the left-hand side of these expressions has been 
suppressed: T v = T v (Q,u). The parameters Tv, Th, etc on the left side of Equations 26 
are characteristics of the scene to be used in integral expressions such as Equation A 6 
(Appendix A) to compute the radiometer output. The dimensions on the left-hand side 
are Kelvin/steradian f 19] which is consistent with this application (i.e. integration by the 
antenna of the radiation per unit solid angle incident from a given direction). 

Notice that the radiometer system bandwidth, Bw, does not appear in these expressions. 
This is reasonable since Equations 26 are parameters of the scene (brightness) and should 
be independent of the measurement device. These expressions are also dimensionally 
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consistent which can be seen by recalling that S pq (Q, u) / r\ has units of brightness and 
that Boltzmann’s constant has units of Joule/Kelvin. 

Writing Equations 26 in terms of electric fields, for example as done in Equations 1 , is 
problematic because of the issues associated with defining Fourier transforms in the case 
of stationary random processes. One can avoid this issue by using the definition of the 
spectral density: 

Spq(u) = F[ Rpq(t) ] (27) 

where F[* \ denotes a Fourier transform and R P q(t) is the correlation function for the 
electric fields incident on the antenna from direction Q: 

Rpq(x) = < e P (t + t) e q *(t) > (28) 

The problem occurs when one tries to rewrite these expressions in terms of the Fourier 
transforms of the fields. In order to accommodate the assumption that the random 
process is stationary, one should use the definition in Appendix C: 

S pq (u,Q) = Lim { <d£p(o,Q)E q *(u,Q)> / 2T} (29) 

T ~} qo 

However, it is appealing to try to preserve the notation in Equations 1 (i.e. the use of 
electric fields on the right-hand side in forms such as < j E v j 2 >). For example, one might 
make a simple substitution of <Ip(Q,n)2q*(Q,u)> for S pq (f2,u) with the understanding 
that the limit in Equation 29 is implied. The difficulty with this approach is that it is 
dimensionally misleading. In particular, < | ‘E v (Q.,v>) j 2 > suggests units of (Volt/meter- 
Hz) 2 whereas S pq (Q,n) in Equation 29 has units of (Volt/meter) 2 /Hz. 

VII. References 


12 



1. J.D. Kraus, Radio Astronomy . McGraw-Hill, 1966. 


a. Chapter 3 

b. Section 6.2. 

2. Ulaby, F.T., R. K. Moore and A.K. Fung, Microwave Remote Sensing . Addison- 
Wesley Publishing Co, 1981. 

a. Section 1-2.2 

b. Section3-2.4 

c. Section 6-7.2 (Equation 6.62) 

d. Section 4-2.1 

e. Section 4-2.2 

f. Section 4-5.3. 

3. Matzler, C., Thermal Microwave Radiation: Applications for Remote Sensing . JET 
Electromagnetic Wave Series, London, UK, 2006 (Section 4.4.3). 

4. Randa. J.. Recommended Terminology for Microwave Radiometry . NIST TN 1551, 
August, 2008; Avail US Government Printing Office; Also see: 
http://wvvw.boulder.nist.gov/div8 1 8/8 1 80 1/Noise/publications/noise pubs.html. 

5. GRSS, http://www.ieee-grss.org/ see “Technical Resources” and then “Radiometer 
Dictionary”, 2008. 

6. Claassen, J.P. and A. K. Fung, “The Recovery of Polarized Apparent Temperature 
Distributions of Flat Scenes from Antenna Temperature Measurements”, IEEE Trans. 
Antennas & Propagation, Vol AP-22 (no. 3), pp. 433-442, May, 1974. 

7. Johnson, J., J.A. Kong, R.T. Shin, D. H. Staelin, K. O’Neill and A.W. Lohanick, 
“Third Stokes Parameter Emission from a Periodic Water Surface”, IEEE Trans 
Geoscience and Remote Sensing, Vol 31 (#5), pp 1066-1080, September, 1993. 

8. Colliander, A., J. Kettunen and M.T. Hallikainen, “Calibration of End-to-End Phase 
Imbalance of Polarimetric Radiometers”, IEEE Trans Geoscience and Remote Sensing, 
Vol 44 (#10), pp 2635-2641, October, 2006. 

9. Skou. N and D. M. Le Vine, Microwave Radiometer Systems: Design and Analysis . 
Artech House, 2006, Chapter 7. 

10. Thompson, A.R., J.M. Moran and G.W. Swenson, lnteferometry and Synthesis in 
Radio Astronomy , New York, Wiley, 1986 


13 



11. Martin-Niera, M., Ribo, S., and Matin-Polegre, A.J., “Polarimetric Mode of MIRAS”, 
IEEE Trans. Geosci. & Remote Sensi., Vol 40, No 8, pp 1755- 1768, August, 2002. 

12. Corbella, I., N. Duffo, M. Vall-llossera, A. Camps, F. Torres, The Visibility Function 
in Interferometric Aperture Synthesis Radiometry, IEEE Transactions Geoscience and 
Remote Sensing, Vol 42 (8), pp. 1677-1682, August 2004, Digital Object Identifier 

10.1 109ATGRS.2004.830641, 

13. Collin, R.E. and F.J. Zucker, Antenna Theory , Vol 1, Chapter 4, McGraw-Hill, 1969. 

14. Rytov, S, Y. Kratsov and V. Tatarskii, Principles of Statistical Radiophysics, 

Elements of Random Fields . Vol 3, Chapter 3, Springer- Verlag, 1989. 

15. Janssen, M.A., Atmospheric Remote Sensing by Mirowave Radiometry , John Wiley 
and Sons, 1993, Section 1.3.1. 

16. Stelzried, C. T., “Microwave Thermal Noise Standards”, IEEE Trans. Microwave 
Theory and Techniques, Vol MTT-16 (no. 9), pp 646-665, 1968. 

17. Born, M and E. Wolf, Principles of Optics . Pergamon Press, 1959. 

a. Section 10.8 

b. Sectionl0.2. 

18. Tsang, L., Xu. P. and Chen, K.S. " Third and Fourth Stokes Parameters in 
Polarimetric Passive Microwave Remote Sensing of Rough Surfaces over Layered 
Media", Microwave and Optical Technology Letters, vol. 50, no. 12, pp. 3063-3069, Dec. 
2008. 

19. Ruf, C.S., C.T. Swift, A. B. Tanner and D. M. Le Vine, “Interferometric Synthetic 
Aperture Microwave Radiometry for the Remote Sensing of the Earth”, IEEE Trans. 
Geoscience and Remote Sensing, Vol 26, No 5, pp597-61 1, September, 1988. 

20. Jordan, E.C. and K.C. Balmain, Electromagnetic Waves and Radiating Systems , 
Section 13.07, Prentice-Hall, 2 ed, 1968. 

21. Jones, D.S., Methods in Electromagnetic Wave Propagation, Oxford Engineering 
Science Series, Clarendon Press, Oxford, Section 9.8, 1979. 

VIII. Appendix A: Traditional Radiometer Response 


14 



The power, W p , available at the output of an antenna with polarization p = {v,h} when 
looking at an unpolarized source with brightness, B(£2,v), is [la, 2e]: 

W p = (1/2) A e If B(Q, v) P np (0 v,) d Q d v (Al) 

where A e is the effective area of the antenna and P n is the normalized antenna pattern. 


Assume that the Rayleigh- Jeans approximation applies: 

B(H, v) = 2 k T (Q, v) / X 2 W/m 2 - sr-Hz (A2) 


and use 

W p = kT A pp Bw (A3) 

on the left side of Equation Al to express the output power in bandwidth Bw in terms of 
the change in temperature, T A PP , of a resistor at the output. Then, Equation A2 can be 
written: 


Bw = A t 


pp DW 


II 


T(Q,v) V 2 P np (Q) dfldv 


(A4) 


If T(Q, v) and P np are independent of frequency over the passband of the measurement 
system and the bandwidth, B w , is narrow compared to the center frequency (usually a 
good approximation for remote sensing observations of thermal sources), then Bw will 
cancel on both sides. Finally, using [lb, 2b]: 


X 2 /Q b 


(A5) 


where £2 B is the antenna solid angle, Equation A4 becomes: 

T A pp = (1/ 0 B ) Jj T(Q) P np (a) df2 (A6) 

Equation A6 is the classic result relating the antenna temperature, T A PP at polarization, p, 
to the source brightness temperature T(Q) [2fj. It applies to a single polarization and 
assumes that the source and antenna parameters that are essentially constant over the 
bandwidth of the measurements. 


IX. Appendix B: Brightness and Electromagnetic Fields 

The objective of this appendix is to establish the relationship between the source 
brightness and the electromagnetic fields incident at the radiometer antenna. 
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A. Time Average Poynting Vector 

The time average of the power per unit area arriving from a direction n (i.e. crossing a 
plane with unit normal, n, located at the antenna) can be expressed in terms of the 
incident electromagnetic fields using Poynting’ s theorem: 

W/A = <E(r,t)xH(r,t)>i • n (Bl) 

where W is the power passing through a plane with area A and unit normal n; and where, 
E and H are the electric and magnetic field intensities, respectively, and <*> t denotes a 
time average, and “x” denotes the “vector” or “cross product”. 

In order to evaluate Equation B 1 , it is convenient to represent the field by their Fourier 
transform and do the analysis on the transforms: 

00 

E(r,t) = J E( r,v) exp(-j2jrvt) dv (B2) 

-00 

To express the fields in the frequency domain we assume that the sources are far away 
and adopt the radio astronomy model by considering the radiation from sources in a small 
solid angle dQ in the direction Q as illustrated in Figure 2. Consider a sphere of radius R 
surrounding the observer (antenna) as shown in Figure 2. The radius R is assumed to be 
large compared to wavelength but small enough that all sources are outside the sphere. 



Figure 2: Geometry for computing incident electromagnetic radiation 
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The electromagnetic fields at the observer can be written in terms of the electric field in 

the aperture dZ where the cone df2 crosses the sphere. Since R is large, one obtains 

approximately [20]: 

d£(r,v) ~ e £(K,v) exp(j kR) / R dZ (B3a) 

d//(r,v) a K x e ®(K,v) exp(j kR) / n R dZ (B3b) 

where k = 2jtu/c is the wave number and K is a unit vector in the direction of propagation 
from dZ to the observation point and e is an arbitrary, unit amplitude, polarization vector. 

To compute the power available at the observer (antenna), the fields from all apertures 
are added and the time average of the Poynting vector of the total field is computed. 
Using the fact that E(r,t) is real (the Fourier transform is hermitean), one can write this in 
the following form: 

oc oc 

<E(r,t)xH(r,t)> t <U/tf)If II p(K,K) Ifk) £*(k’) exp[j 2ji(v-v’)t] ... 

-00 -00 

. . . exp[j27i(k-k , )R] dv dv’ dZ(K)dZ(K’)/R 2 > t (B4) 

where 

p(K,K’) = [ e k x (K’ x e k 0 ] (B5) 

The fields propagating in the direction K in each aperture, dZ, in Figure 2 are composites 
of radiation from many sources (all the sources in the solid angle dO) and will be 
modeled as a random process. We consider the case of “natural” or “thermal” radiation 
in which case it will be assumed that the radiation from different apertures (directions, K) 
is uncorrelated and that the process is ergodic and stationary [14]. 

Ergodicity permits replacement of the time average in Equation B4 with the statistical 
(ensemble) average. Because the process is stationary and the fields in different patches 
are uncorrelated, one can express the average as: 
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<'£(K, u) £*( K\ d’)> = 


Se(K, u) 5(v - v’) 6(K - K’) 


(B6) 


and Equation B4 becomes: 

<E(r,t)xH(r,t)> = (i in)\ I K Se(K, V) dv dI(K)/R 2 (B7) 

Recognizing that dI(K')/R 2 = dQ where Q denotes solid angle, one can write: 

<E(r,t)xH(r,t)> = (1 /tJ J K(Q) S e (H, u) dv dD (B8) 

The function, S E (fi, u), in Equation B8 is the power spectral density of the electric fields 
arriving from a given direction (0,(p) per unit solid angle, fl It is not the power spectral 
density of the electric field at the antenna but, rather, a “directional” power spectral 
density. See Appendix C for a discussion of the spectral density in this context. 

Also, note that Equation B8 is independent of distance, R, the radius of the sphere in 
Figure 2. This is to be expected because the fields in Equations B3a-b are radiation fields 
and the radiated power is conserved. 

B. Brightness and Electric Field 

From Equation B8, the dimensions of S E (fi, u)/q are seen to be “power/(unit area - 
steradian - Hz)” which are the same as those of brightness, B(f2, u) in Appendix A. In 
addition, comparison of Equations A1 and B8 shows that S E (D, u)/q and B(Q, u) play the 
same role in describing the incident power. The difference is that the antemia has not 
been included in Equation B8. Doing so would result in identical expressions and allow 
us to make the identification: 

B(Q, i>) = W/(m 2 -Hz-sr) 

= S e (A u)/r| (B9) 

Notice that S E (0, u)/q represents the total power available (i.e. from both polarizations) 
and in this sense differs from S pq (D, n)/q defined in the text. The advantage of the 
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approach adopted in the text is that it handles polarization explicitly (see Section V). In 
the case of unpolarized radiation: 

S E (0,n) = S pp (Q,u) + S qq (Q,n) (BIO) 

X. Appendix C: Fourier Transforms 

The transition from Equations 6 to 8 is obtained by recognizing that the given 
assumptions (incident fields that are stationary and uncorrelated in direction) imply that 
the expected value of the integrals over frequency and direction is equivalent to the sum 
over only those terms for which v = v’ and Q = IT and zero otherwise. This result is 
formally equivalent to using Equation 7 in Equation 6. 

The function, S pq (o, Q ), on the right hand side of Equation 7 is the cross-spectral density 
in the classical sense (Fourier transform of the cross-correlation function): 

Spq(u) = F 1 Rpq(r) ] = F ( <e p (t + x) e q *(t)>l (Cl) 

where the dependence on Q has been suppressed. One could try to write Equation Cl in 
terms of the Fourier transforms of e p (t) and e q (t) in order to relate it to the incident 
electromagnetic fields. However, the assumption of stationarity presents some 
mathematical problems for the existence of Fourier transforms (because the functions 
must extend to infinity). 

This issue is formally handled by defining the cross-spectral density S pq (u,n) in a 
limiting sense using the Fourier transform of functions with a defined time window in the 
limit as the width of the window, T, becomes infinite [17b, 21]: 

5 pq (u,Q) = Lim {<£p(n,f2)£' q *(u.n) >/2T } (C2) 

T -> ao 
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where bold face, £ p (u, Q ), indicates Fourier transform of e p (t) defined over the finite 
interval, -T < t < T. Asa check, the dimensions on the right-hand side of Equations C2 
and C2 are equal and that the dimensions of S pq (o, £2 )/ q are those of “brightness”: 
W/(m 2 -Hz-sr). 
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